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We consider a strongly interacting one-dimensional (ID) Bose-Fermi mixture confined in a har- 
monic trap. It consists of a Tonks-Girardeau (TG) gas (ID Bose gas with repulsive hard-core in- 
teractions) and of a non-interacting Fermi gas (ID spin-aligned Fermi gas), both species interacting 
through hard-core repulsive interactions. Using a generalized Bose-Fermi mapping, we determine the 
one-body density matrices, exact particle density profiles, momentum distributions and behaviour of 
the mixture under ID expansion when opening the trap. In real space, bosons and fermions do not 
display any phase separation: the respective density profiles extend over the same region and they 
both present a number of peaks equal to the total number of particles in the trap. In momentum 
space the bosonic component has the typical narrow TG profile, while the fermionic component 
shows a broad distribution with fermionic oscillations at small momenta. Due to the large boson- 
fermion repulsive interactions, both the bosonic and the fermionic momentum distributions decay 
as Cp~* at large momenta, like in the case of a pure bosonic TG gas. The coefficient C is related to 
the two-body density matrix and to the bosonic concentration in the mixture. When opening the 
trap, both momentum distributions "fermionize" under expansion and turn into that of a Fermi gas 
with a particle number equal to the total number of particles in the mixture. 



PACS numbers: 05.30.-d,67.85.-d,67.85.Pq 
I. INTRODUCTION 



Dilute degenerate Bose-Fermi (BF) or Fermi-Fermi 
mixtures have been realized over the past few years in 
several experiments by trapping and cooliiig either gases 
made of mixed alkali-atom isotopes [l|, [2i H 0, S Q i or 
of imbalanced two-spin fermionic atoms (7|, |8|. In the 
latter, a superfluid paired core is surrounded by a shell 
of normal unpaired fermions. For BF mixtures a mean- 
field approach predicts that the boson-fermion coupling 
can lead to quantum phase transitions and, in partic- 
ular, that boson-fermion repulsion can induce a spatial 
demixing of the bosonic and fermionic components when 
the interaction energy overcomes the kinetic and confine- 
ment energies @. The instability leading to this phase 
separation has been studied in all dimensions both for 
homo gen eous 0, and for harmonically-trapped mix- 
tures [111, [13, [Q ■ As an illustrative example, we show 
in Fig. [T] the bosonic and fermionic density profiles 
nsix) and np^x) as obtained from a two-fluid mean- 
field model when a BF mixture with equal masses m 
is confined in a ID harmonic trap and phase separation 
occurs as the Pauli pressure overcomes the boson-boson 
repulsion [l3|. In the model, the particles experience con- 
tact boson-boson vbb{x) — gsB^ix) and boson-fermion 
vbf{x) = gBF^ix) interactions. The quasi-lD interac- 
tion strengths gsB and gsF, which can be expressed in 



terms of the 3D scattering lengths [IJI , define in turn the 
adimensional coupling constants jbb — fi^gsB / {mriB) 
and jBF = ^gBF I {fnnp). These parameters measure 
the ratio of the interaction to kinetic energies. Rather 
counter-intuitively, in ID, the weakly interacting regime 
corresponds to the large-density regime. The mean-field 
description applies well to local observables in the weakly 
interacting regime, however, because phase fluctuations 
do play a major role in ID, the off-diagonal long-range 
order is lost and the one-body density matrix decays as 
a power-law, even at zero temperature, with an expo- 
nent governed by the strength of the interactions (see 
e.g. [il). 

In this work, we focus on the strongly interacting 
regime ^bb 3> 1, 7b_f 3> 1 which is beyond the regime 
of validity of the mean- field treatment. This strongly in- 
teracting regime is not far from experimental access. In- 
deed the increasing sophistication of experimental tech- 
niques allows to realize traps so tight that the atomic 
dynamics is essentially one-dimensional (Tg} and to drive 
the effective ID coupling strength to very large val- 
ues [itI . [ist . Previous theoretical studies of strongly 
interacting BF mixtures include the extension of the 
Bethe Ansatz method developed by Lieb and Lininger for 
bosons to the homogeneous BF mixture (sec e.g., [19] and 
references therein), methods from conformal quantum 
field theory |20j, a numerical analysis of a harmonically 
trapped Bose-Bose mixture using a multi-configuration 



2 



0.016 ^ ^ J- ^ ^ 

0.012 - 

^ o.oos - 

0.004 - - 

n I — ^ ' LJ ' ^ — 

-400 -200 200 400 

X 

FIG. 1: (Color online) Bosonic (red solid line) and fermionic 
(x20, green dashed line) density profiles (in units of l/£o and 
normalized to the total number A'^ of particles) as a func- 
tion of the distance x from the trap center (in units of ^o) 
as predicted by a two-fluid mean-field model [13 ] for a BF 
mixture confined in a ID harmonic trap. The length scale 
unit £o is the harmonic oscillator length as defined in Sec. 
im The boson-boson and boson-fermion interaction strengths 
have been choosen such that jbb ~ 1.5x 10~* and 7sf = 0.63 
(see the text for the definition of "/bb and "/bf)- 



the hard-core boson-fermion repulsion. We also study the 
ID expansion of the mixture when the trapping external 
potential is turned off and the BF mixture is released 
from the trap and is allowed to expand in a ID waveg- 
uide. The net result of this expansion is to "fermionize" 
the initial fermion and boson momentum distributions: 
once the expansion is completed, the shape of the re- 
sulting distributions in momentum space is the same as 
the shape of the initial ones in real space, and they both 
present a number of peaks equal to the total number of 
particles. Our work thus generalizes the results obtained 
for a pure TG gas in [2^, [2^. 

The paper is organized as follows. Section |TT] intro- 
duces the model for the BF mixture under consideration. 
In Sec. mil we derive the bosonic and fermionic one-body 
density matrices and, in Sees. IIVI and fVl respectively, we 
compute analytically the corresponding spatial densities 
and momentum distributions. The effects of the expan- 
sion dynamics on the momentum distributions when the 
trap is turned off are described in Sec. IVII Finally Sec- 
tion lVIII offers some concluding remarks and perspectives. 



time-dependent Hartree method [2l| and the exact so- 
lution for the many-body wavefunction using a general- 
ized Bose-Fermi mapping method valid when the Tonks- 
Girardeau (TG) limit jbb oo and 7^^? — )■ oo is ob- 
tained for both species [2^ . 

In this paper we analyze several equilibrium and dy- 
namical properties of the inhomogeneous Bose-Fermi 
mixture in the TG regime. We use the TG many-body 
wavefunction determined in Ref. [2^ to obtain the exact 
density profiles, momentum distributions and behaviour 
under longitudinal expansion of a mixture of bosons and 
fermions (with equal masses m) trapped by the same ex- 
ternal harmonic potential. The bosonic component of 
the mixture is made of identical impenetrable bosons 
(TG gas) and the fermionic component is made of non- 
interacting identical spin-polarized fermions. Finally the 
boson-fermion interaction is characterized by a point-like 
infinite hard-core repulsion. As a whole the BF mix- 
ture under consideration is thus made of impenetrable 
particles. In the TG limit, at variance with the mean- 
field predictions, this system does not display any phase 
separation in real space: the boson and fermion density 
profiles are proportional to each other, they extend over 
the same region of space and both present a number of 
density peaks equal to the total number of particles in 
the trap. The bosons and fermions display however dif- 
ferences in their momentum distributions, especially at 
small momenta where the bosonic momentum distribu- 
tion is sharply peaked around p ^ while the fermionic 
one, refiecting the Fermi sphere occupancy, is broad. 
Quite remarkably, we also find that both momentum dis- 
tributions have the same p~'^ asymptotic behaviour as 
the pure hard-core Bose gas [H, . We relate the co- 
efficient of this power-law decay to the 2-body density 
matrix of the system. For the fermionic component, this 
" bosonization" of the momentum distribution is due to 



II. THE MODEL 

We consider a mixture of Nb identical bosons with TG 
interactions and Np identical non-interacting fermions 
having equal masses m and trapped in the same ID 
harmonic potential with trapping frequency wq. Every 
boson-fermion pair is subjected to TG (repulsive) inter- 
actions. Throughout the paper we will use dimension- 
less variables, expressing all spatial variables in units of 
the harmonic oscillator length £0 = ^Jfi/muio, all mo- 
menta in units of the harmonic oscillator momentum 
Po — ^/^o = VmhujQ, all energies in units of the har- 
monic oscillator energy fujQ and time in units of uJq^. 
By convention, we will collectively denote the space co- 
ordinates by X — {xi,...,xn) understanding that Xi is 
a scaled bosonic variable provided i G B = {1, Nb} 
while it is a scaled fermionic one provided i ^ T — 
{Nb + l,...,iV_B + A^F = N}. The same convention 
will apply for scaled momenta denoted collectively as 

P = (pi, ...,PAr). 

In the Tonks-Girardeau limit, when ^ 00 and 
^BF — *■ 00, the interactions are appropriately incorpo- 
rated in the model by imposing boundary constraints to 
the total wavefunction. In this case of infinite repulsive 
contact interactions, and in terms of the scaled variables, 
the total Hamilton differential operator is then a sum of 
one-body Hamilton differential operators H = X^iLi ^(*) 
where h(i) = {—d^. + xf)/2 in space representation. We 
now demand that '^{X) = whenever Xi = Xj for i ^ j. 
Using the Fermi-Bose mapping theorem [l^], the total 
wavefunction satisfying these constraints is constructed 
from a Slater determinant and reads 

^{X) = A{X)-^d{X). (1) 
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Here the antisymmetrizor is given by 

MX) = [n sgn(a;, - x^)] [ n n «gn(a 



(2) 



and the determinantal wave function is 



N 



(3) 



where P runs over aU A^! possible permutations, including 
boson-fermion exchanges, e(-P) = ±1 is the sign of the 
permutation and the {ui\ are the single-particle orbitals 
for the given external potential. We recall that the choice 
^ for the mapping function A{X) is not unique as in the 
limit jBB oo and ^bf oo the ground state has a 
large degeneracy. We choose here the solution which has 
the same cusps as the lowest-energy solution for finite 
interaction strength [2^ ]. 

For the harmonic potential the groundstate wave func- 
tion is thus obtained by taking the rescaled Hermite- 
Gaussian orbitals = {^/^^ 2"n\)~^/'^ ^^Hn{x). 

As already shown in Ref. 0|, the full many-body 
groundstate wavefunction can then be written as 



fo{X) 



with 



i>3 



and the iV-dependcnt 



(4) 

j-xi\ 

(5) 

cjv = 



2^(^-1)74 



x\X\lzl 



normalization constant 

1/2 

. One can easily check 



that indeed fulfills the constraints imposed by quan- 
tum statistics and the infinite repulsive contact interac- 
tions. It is duly symmetric under the exchange of any 
two bosons, it is duly antisymmetric under the exchange 
of any two fermions and it vanishes as soon as any two 
particles are at the same location in space. Furthermore 
it is straightforward to check that l^'oP is normalized to 
unity, i.e. /dX|*o(^)P = 1- 



III. ONE-BODY DENSITY MATRICES 

The one-body bosonic density matrix p"^ is obtained 
by tracing out all particles except one boson. Since \l/o is 
symmetric under the exchange of bosons, it can always 
be written as 



p%\x,y) = NbJ dX' ^oix,X')%{y,X') 



(6) 



where X' is a shorthand for all other variables Xi except 
the first bosonic one (i.e. i = 2, . . . , N). Since is nor- 
malized to unity, p'^^ is normalized to the total number 



of bosons Nb- It can be shown that p^-* is in fact propor- 
tional to the one-body density matrix of a pure TG gas 
of TV bosons namely p%\x,y) = {NB/N)pTG{x,y) 



|27| so that we finally get 

PB^x^y) 



2 -(x=+j.^)/2dgtr(x,y) (7) 



TTiV! 



where T{x,y) is a square matrix of order (N — 1) with 
entries 



^jkix,y) 



^r(j)r(fc) 



dtf 



x^t\\y-t\t^+'^-\ 
(8) 

and r(j) = (j — 1)1 is the Gamma function. 

The determination of the fermionic one-body density 
matrix {x, y) follows the same rationale, namely inte- 
grate over all particles in the mixture except one fermion 
(which can be chosen to be the last one, i.e. j = N). It 
has been calculated in [l^, in the case of a homoge- 
neous mixture and the calculation proves more complex. 
We have extended the result derived in [2^ to the har- 
monically trapped mixture and we find: 

P^^^(-.^)-4SSr^e-(^^+^^)/^A(x,,) (9) 



/i^m {N - 1)! 



where 



^{x, y)^J ^e-'^^- det [r{x, y)e'^ + A(x, y)] 

(10) 

with A(x, y) being a square matrix of order (N — 1) with 
entries given by the Gaussian integrals 



^3k{x,y) = 



nr{j)T{k)J_ 



die"* {x-t){y-t)P 



+k-2 



(11) 

The integration over the phase (j) in (fTO|) ensures that only 
the terms involving Nb bosons are picked up and the fac- 
torials in ([9]) take care of multiple-counting. Equation ([9]) 
could not be reduced further to a simpler analytical for- 
mula but the preceding expressions allow for an efficient 
numerical computation of p^p' . 

It is possible to study the large-distance off-diagonal 
behaviour of the one-body density matrix by factoring 
out the Gaussian functions imposed by the trap. In prac- 
tice, we rescale the bosonic and fermionic one-body den- 
sity matrices by the square root of their corresponding 
density profiles, i.e. we define the one-body correlators 
as in l28ll. 



9^'\x,y) 



jP^x^y)_ 
^/n{x)n{y) 



(12) 



In Fig. [21 we display the off-diagonal behaviour of 

gg\x,y) and gp\x,y) as obtained for y — ~x. Fig. [2la, 
shows the results for a Bose-Fermi mixture with = 20 
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particles when the boson concentration is varied while 
Fig. [2t and Fig. [2Ji show the same quantity for a pure 
TG gas and a pure Fermi gas. At large distances, the 
fermionic one-body correlator of the BF mixture displays 
a power-law decay |a;|~" modulated by typical fermionic 
oscillations. The wave vector of the oscillations is close to 
(but not exactly equal to) 2kp = 27rnF(0) as evidenced 
in Fig. The power-law exponent a depends on the 
boson concentration in the mixture as is illustrated in 
Fig. d). For Nb = Np ^ N/2 the exponent is a ~ 2.6, 
i.e. larger than that of a pure Fermi gas, where a — 1 
(for a homogeneous gas in the thermodynamic limit), as 
well as that of a pure TG gas, where a = 1/2 (again 
for a homogeneous gas in the thermodynamic limit). We 
therefore find that the Bose-Fermi interactions strongly 
affect the one-body correlations at large distance. 



IV. DENSITY PROFILES 

Since the square of an anti-symmetric operator is the 
unity, the bosonic and fermionic density profiles n b (x) — 

P'b (a;, x) and np^x) — p'p (a;, x) are exactly the same up 
to normalization factors. Since p'b proportional to 
Ptg, we have that nB(x) and npix) are both propor- 
tional to the spatial density nToix) = ptg{x,x) of a 
harmonically trapped TG gas made of = Nb + Np 
bosons [22] : 



a) 



nB{x) np{x) nToix) 



Nb 



Np 



N 



(13) 



The main consequence of this result is that strongly inter- 
acting spin-polarized fermions and impenetrable bosons 
do not display any phase separation (or spatial demix- 
ing). Both the fermionic and the bosonic density profiles 
present a number of peaks equal to the total number of 
atoms and the position of the peaks is exactly the same 
for bosons and fermions. As an illustrative example, we 
plot in Fig.[3]the fermionic density profiles obtained for a 
BF mixture with a fixed total number of A'^ = 7 particles 
when the number of fermions is increased from Np —1 
up to 7. 

We also note in passing that for a purely fermionic en- 
semble, the Hamilton differential operator does not dis- 
tinguish positions and momenta up to scaling. Hence, for 
noninteracting fermions the functional dependence of the 
momentum distribution vs momentum p is exactly the 
same as the functional dependence of the density profile 
vs X (dashed curve in Fig. [3] for the case of 7 fermions). 
In the next section we discuss how this functional depen- 
dence is affected by the interaction with the bosons. 



V. MOMENTUM DISTRIBUTION 

Using the scaled variables, the momentum distribution 
is defined as the Fourier transform of the one-body den- 
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FIG. 2: (Color online) Plot of the off-diagonal behaviour 
of the fermionic and bosonic one-body density correlators 
gp\x, —x) and gg\x, —x), Eq. (|12p . as a function of the 
spatial coordinate x (in units of £o). Panel (a) shows g^p^ 
for the BF mixture (in log-log scale) with A'^ = 20 par- 
ticles and three different bosonic concentrations Nb/N = 
.75 (-I-), .5 (*), .25 (□). The comparison to the fitting function 
~ \x\~°' (continuous, short-dashed, dotted) allows to extract 
the slope. Panel (b) shows the results for the power-law ex- 
ponent a obtained at large distances x for the mixture with 
= 20 particles when the boson concentration TVs /N is var- 
ied. Panel (c) shows the log-log plots of g^ -' for a pure TG 

gas with A'' — Nb ~ 20 particles (*) and of g^^' for the pure 
Fermi gas with A*' — Nf = 20 particles (+). For the pure TG 
gas, the comparison to the fitting function ~ |a^|~" (double- 
dotted line) gives a — 0.44, a value close to the one found for 
the homogeneous case where a — 0.5. For the pure Fermi gas, 
the comparison to the fitting function ~ sm{2kFx) /\x\" (dot- 
dashed line) gives a — 0.93, a value close to the one found for 
the homogeneous case where a — 1. Panel (d) is the same as 
panel (c) but in linear scales. 



sity matrix according to: 

nip) ^ — / dx dyp^^\x,y)e 

J-oo J-oo 



(14) 



(1) 



The case for bosons is easily cleared. Indeed, since 
and ptg are proportional, UBip) is identical to nTcip) 
(up to normalization). But the case for fermions is more 
involved because the expression ^ for the density matrix 
could not be reduced to a simpler analytical formula. 



FIG. 3: (Color online) Fermionic density profiles npix) (in 
units of l/£o) versus x (in units of ^o), as obtained for a BF 
mixture with a fixed total number of = 7 particles when 
the number of fermions is increased from Nf ~ 1 (solid line) 
up to Nf — 7 (double-dashed line). Normalizing each curve 
by Nf/N would collapse them on a same and single curve. 



FIG. 5: (Color online) Fermionic momentum distribution 
^f('p) (in units of 1/po) versus p (in units of po), at fixed 
number of fermions Nf = 4 when bosons are gradually added. 
The curves range from Nb = (solid line) down to Nb = 3 
(double-dashed line). As one can see, adding more and more 
bosons smoothens the fermionic oscillations. 




FIG. 4: (Color online) Fermionic momentum distribution 
nF(p) (in units of 1/po) versus p (in units of po), as obtained 
for a BF mixture with a fixed total number of A = 7 particles 
when the number of fermions is raised from Nf =1 (solid line) 
up to Nf =7 (double-dashed line). As fermions are added, 
Nf "fermionic" oscillations develop in the momentum distri- 
bution. 



and we have resorted to numerical computations. Fig. 
|4] shows our results for a BF mixture with a fixed total 
number of A'^ = 7 particles when the number of bosons is 
increased from Nb — up to 6. Fermionic oscillations are 
suppressed and the tails of the distribution become more 
prominent. The large-p behaviour of the distributions 
will be analyzed here below. 

To evidence how the presence of bosons can influence 
the spread of np{p)^ we have also considered an ensem- 
ble of Np — 4 fermions and we have gradually added 
bosons to it. In Fig.[5]we have plotted the corresponding 
fcrmionic momentum distributions when Nb is increased 
from up to 3. Adding more and more bosons smoothens 
the fermionic oscillations while the tails of the distribu- 
tion get enhanced and the distribution broadens. This 
effect can be understood by looking at the asymptotic 
behaviour of the fermionic momentum distribution when 
p oo. 



To this purpose, we use the alternative expression 

np{p)^Np f dY\i'o{Y,p)\^ (15) 



where y is a shorthand for the collection of all Xi vari- 
ables except the last fermionic one xn (i-e. i — 1, ■ ■ ■ N ~ 
1). Here >I'o(F,p) stands for the Fourier transform, eval- 
uated at p, of the many-body wavefunction ^'o(F, xat) 
with respect to its last fermionic variable xn 



"foiY^p) 



dx 



N 



'2ti 



^o{Y,XN)e- 



(16) 



We now use the mathematical result that for any f{z) 



F{z), where F{z) is a regular function and 



a > 0, a 7^ 1, 3, 5, • • •, we have 



lim 



fL/2 
J -L/2 



:2F(zo)cos(^) F(a). 



_ (17) 

Following the procedure outlined in Ref. [2J|, we next 
obtain the large-p asymptotics of ^o{Y,p) 

*o(r,p) ^ ^£^/^(r)if(r)e-^f- (i8) 



where /q is given by ^ but with the fermionic index now 
only running in J-' = J-\{N} and 



K{Y) = Y.^ 



— ipxj 



n n \Xj - Xl\{Xn- Xj). 



(19) 

The momentum distribution then follows by squaring 
and integrating over Y. After some manipulations, we 
get np{p) ^ Cp~'^ where the constant C is given by 



2NfNb 2 



dxe-^'' D{x). 



(20) 
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The function D{x) reads 



N 

D{x) = /dC/ Jl [e^ 

i=3 



i-1 

i=3 



(21) 



U being the shorthand for the coUection of aU Xi for 
2> < i < N . This expression can be compacted to a de- 
terminant form by noting that D(x) = (iV — 2)! det G{x) 
where G{x) is the square matrix of order [N — 2) with 
entries 



Gjkix) = / dte-' [x-tft^ 



+ fe-2 



(22) 



Each of these matrix elements can be evaluated using 
Gamma functions [131 ■ As a net result the coefficient C 
can alternatively be written as 



C 



2NfN, 



B 



(TV -2)! 4 



N 



dx e 



detG{x). (23) 



The asymptotic result n{p) ^ Gp^"^ captures well the 
large-p behaviour of np{p). This is illustrated in a log- 
log graph in Fig. [SJ where a good agreement is displayed 
between the asymptotic result and the full calculation of 
the momentum distribution p4p . 




FIG. 6: (Color online) Momentum distribution npip) (in 
units of 1/po) as a function of p (in units of po) for the case of 
a BF mixture made of Nf = 2 fermions and Nb ~ 5 bosons. 
The plot shows the comparison between the analytical asymp- 
totic prediction Cp~*, with C being given by (|23p . (dashed 
line) and the exact result computed from Eq. (|14p (squares). 

The asymptotic behaviour of the fermionic mo- 
mentum distribution in the BF mixture under considera- 
tion is the same as that of a TG gas [l^ . It is straightfor- 
ward to show that (|20| or (p3)) can indeed be rewritten 
as 



C = lim p'^np{p) — 



NFNf_ 



N{N - 1) p->oo 



lim pWcip) (24) 



where nTcip) the momentum distribution for a har- 
monically trapped TG gas of N identical bosons. The 
factor Np /N ensures that np [p) is indeed normalized to 
the total number Np of fermions in the system, while 



the remaining factor Nb / {N — 1) gives the weight of the 
bosonic contribution at large p. We can also notice that, 
for a BF mixture with a fixed number N of particles, the 
effect of the bosons on the fermionic distribution will be 
maximized when the number of bosons equals the num- 
ber of fermions, i.e. Np = Nb- 

Interestingly enough, the coefficient G can be related to 
the two-body density matrix p'^' . The two-body density 
matrix is the same for both the bosonic and fermionic 
components as it does not depend on the sign of the 
many-body wavefunction, and coincides with the two- 
body density matrix of an ideal Fermi gas in the same 
external potential. It is defined as 

(x, V- x\ v')=N{N - l)[ dZ^oiZ, x, y)%{Z, x' , y'), 



(25) 

where Z is, a, shorthand for the collection of all Xj except 
the last two fermionic variables (i.e. j = 1, • • • , — 2). 
Using the explicit definition of p*-^) for the mixture we 
can rewrite (PU)) as 



C 



2NpNb 



nN{N- 1) 



a. ,i„ fi^i^^ (26, 

^ x\x"^x \X — X \\X — X \ 



Using (|24p . we also obtain the large-momentum be- 
haviour of a pure TG gas as 



hm p riTGip) = - 

p — ^oo TT 



da; lim 



p^^'>{x\x;x",x) 



*x \X — X \\X 



(27) 

This expression is the TG limit extension of the results 
found by Olshanii [13] in the case of a homogeneous gas 
of hard-core bosons on a wire of length L with finite 
interactions described by the ID s-wave scattering length 
a\p)i namely 



p('HO,0;0,0) 



{plhY 



(28) 



'■YD 



Here p is given in full units and not rescaled to the har- 
monic oscillator ones. Furthermore, the distribution of 
discrete momenta rip is normalized as X)j^-oo "■27rSj/L = 
N . In the TG limit, both the ID scattering length and 
the diagonal element of the two-body correlation func- 
tion p(^)(0, 0; 0, 0) vanish; our expression (P7|) gives the 
corresponding limiting value. 



VI. EXPANSION 

We now turn to the dynamical evolution of the har- 
monically trapped BF mixture when the trap is switched 
off by turning down to zero the trap frequency according 
to Lo(t) — uJof{t) with some known function f(t). In a 
real 3D experiment, a ID expansion could be generated 
by turning off only the longitudinal confinement and the 
mixture would then expand in a ID geometry. This dy- 
namics can be described exactly by using time-dependent 
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coherent states [29| in close analogy with the dynamics 
of a TG gas [l^ . To this purpose we introduce a scahng 
transformation, acting on both the spatial and temporal 
coordinates, which provides an exact solution describing 
this expansion in terms of the time-dependent orbitals 
Uj{x, t) with energy Ej 



Uj(x, t) 



1 



(29) 



The scaling factor hit) obeys the ordinary differential 
equation h + ph = b^^ with initial conditions 6(0) = 1 
and 6(0) = 0. Finally the rescaled time r is given by 
r(t)=/>'[6(f)]-2.^ 

By exploiting the time-dependent Bose- Fermi mapping 
[30| , we construct the many-body wave-function in terms 
of the orbitals Uj{x,t), and hence we obtain the time 
evolution of the one-body density matrix in the following 
scaling form (26l | 



p«(x,2/;t) 



pW(x^;0) 

m 



exp 



ib{x'^ — y^) 
26(1) 



(30) 



The phase acquired by the one-body density matrix is 
crucial for the time evolution of the momentum distribu- 
tion, which is determined by 



nip,t) = 6/d2;dyp(i)(a;,y;0)e-'''''(^'-«')/2-M^- 



v) 



(31) 

We now analyze the case of a sudden turn-off of the 
trap, i.e. f{t) = 1 — 8(t), Q being the step function. The 
scaling parameter b{t) is then found to be b{t) — Vl + t^. 
In Fig. [7] we compare the time evolution of the momen- 
tum distributions obtained for Nb = 7 bosons in the 
absence of fermions, for Np = 2 fermions mixed with 
Nb — 5 bosons and for Np — 7 fermions in the absence of 
bosons. The agreement among all three momentum dis- 
tributions at large times (bottom-right frame) provides 
clear evidence for the dynamical fermionization at work 
in the BF mixture, hence generalizing the results for a 
pure TG gas obtained in [2^ . 



VII. SUMMARY AND CONCLUDING 
REMARKS 

We have determined exactly the groundstate proper- 
ties of a strongly interacting mixture composed of non- 
interacting fermions and hard-core point bosons with mu- 
tual repulsive point-like hard-core interactions, subjected 
to a ID harmonic confinement. Using a generalized Bose- 
Fermi mapping (2^ designed to extend the Bethe Ansatz 
approach [l9| to inhomogeneous confinements, we have 
explicitly determined the bosonic and fermionic one-body 
density matrices. The bosonic one-body density matrix 
is proportional to that of a pure Tonks- Girardeau gas 
whereas the fermionic one can be given a form suitable for 




-6 -4 -2 2 4 e 



FIG. 7: (Color online) Scaled momentum distribution n{p) 
(in units of 1/po) versus p (in units of po) of an expanding 
mixture at different times (in units of ojq^) as indicated in 
the panels. The different lines correspond to Nb = 7 bosons 
(continuous red lines), Np = 2 fermions mixed with Nb = 5 
bosons (dashed green lines), and Np — 7 fermions (double- 
dashed black lines). The distributions have been rescaled by 
the total number A'' of particles in order to be normalized to 
unity. 



practical calculations by expressing it in terms of a one- 
dimensional integral involving the determinant of known 
special functions. Knowing the one-body density ma- 
trices, we have studied in detail the equilibrium density 
profiles, the momentum distributions and the behaviour 
of the latter under one-dimensional expansion. Concern- 
ing the density profiles, at variance with the mean- field 
predictions, no phase separation occurs in the strongly in- 
teracting regime. Both the fermionic and bosonic density 
profiles are proportional to each other and display a shell 
structure with oscillations where the number of peaks is 
equal to the total number of particles in the mixture. 
While the bosonic momentum distribution keeps propor- 
tional to that of a pure TG gas, we have found that the 
fermionic momentum distribution shows instead features 
of a Fermi sphere at small momenta and slowly decaying 
tails, like Cp~*, at large momenta. This is a character- 
istic feature due to the interactions of the fermions with 
their bosonic partners. We have also determined analyt- 
ically the coefficient C in terms of the two-body density 
matrix. This coefficient might be measured in an expan- 
sion experiment as suggested by Tan ^31] in the 3D case. 
We have also studied the behaviour of the mixture after 
a sudden turn-off of the harmonic trap. We have found 
that the momentum distribution of the mixture "fermion- 
izes" during the expansion: its long-time limit shape is 
that of a non-interacting Fermi gas with a total number 
N — Nb + Np of particles. Our predictions could be 
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verified by the on-going experiments on ultracold atomic 
mixtures loaded in optical lattices in the limit of low fill- 
ing factor. 
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